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Abstract 

Second- and third-order results are presented for the structure functions of charged-current deep- 
inelastic scattering in the framework of massless perturbative QCD. We write down the two-loop 
differences between the corresponding crossing-even and -odd coefficient functions, including 
those for the longitudinal structure function not covered in the literature so far. At three loops 
we compute the lowest five moments of these differences for all three structure functions and 
provide approximate expressions in Bjorken-x space. Also calculated is the related third-order 
coefficient-function correction to the Gottfried sum rule. We confirm the conjectured suppression 
of these quantities if the number of colours is large. Finally we derive the second- and third-order 
QCD contributions to the Paschos-Wolfenstein ratio used for the determination of the weak mixing 
angle from neutrino-nucleon deep-inelastic scattering. These contributions are found to be small. 



1 Introduction 



Structure functions in deep-inelastic scattering (DIS) are among the most extensively measured 
observables. Today the combined data from fixed-target experiments and the HERA collider spans 
about four orders of magnitude in both Bjorken-jc and the scale Q 2 = —q 2 given by the momentum 
q of the exchanged electroweak gauge boson [1]. In this article we focus on the ^-exchange 
charged-current (CC) case, see Refs. [2-4] and [5-8] for recent measurements in neutrino DIS and 
at HERA. With six structure functions, F 2 W , F^ and F^ , this case has a far richer structure 
than, for example, electromagnetic DIS with only two independent observables, F2 and F L . 

More detailed measurements are required to fully exploit the resulting potential, for instance 
at a future neutrino factory, see Ref. [9], and the LHeC, the proposed high-luminosity electron- 
proton collider at the LHC [10]. Already now, however, charged-current DIS provides important 
information on the parton structure of the proton, e.g., its flavour decomposition and the valence- 
quark distributions. Moreover, present results are also sensitive to electroweak parameters of the 
Standard Model such as sin Oiv, see Ref. [11], and the space-like W-boson propagator [12]. As 
discussed, for example, in Refs. [13-16], a reliable determination of sin 0^ from neutrino DIS 
requires a detailed understanding of non-perturbative and perturbative QCD effects. 

The perturbative calculations for the unpolarised structure functions in DIS have almost been 
completed to the next-to-next-to-leading order (NNLO) of massless QCD. These results include the 
splitting functions, controlling the scale evolution of the parton distributions, to the third order in 
the strong coupling constant oc s [17, 18], as well as the hard- scattering coefficient functions for F\, 
F2 and F3 to second order in a s [19-23]. For the longitudinal structure function Fl = F2 — 2xF\ the 
third-order coefficient functions are required at NNLO. So far these quantities have been computed 
only for electromagnetic (photon-exchange) DIS [24,25]. In fact, it appears that even the second- 
order coefficient functions for the charged-current Fl have not been fully presented in the literature. 

It is convenient to consider linear combinations of the charged-current structure functions F fl M ' ± 
with simple properties under crossing, such as Fa P±Vp (a = 2, 3, L) for neutrino DIS. For all 
these combinations either the even or odd moments can be calculated in Mellin-/V space in the 
framework of the operator product expansion (OPE), see Ref. [26]. The results for the third-order 
coefficient functions for the even-TV combinations F 2 [ +yp can be taken over from electromagnetic 

DIS [24, 25]. Also the coefficient function for the odd-TV based quantity F^ p+yp is completely 
known at three-loop accuracy, with the results only published via compact parametrizations so 
far [27]. For the remaining combinations F 2 P Vp and F^ p Vp , on the other hand, only the first five 
odd and even integer moments of the respective coefficient functions have been calculated to third 
order in Ref. [28] following the approach of Refs. [29-31] based on the Mincer program [32,33]. 

The complete results of Refs. [24, 25, 27] fix all even and odd moments TV. Hence already 
the present knowledge is sufficient to determine also the lowest five moments of the differences 
of corresponding even-TV and odd-TV coefficient functions and to address a theoretical conjecture 
[34] for these quantities. Furthermore these moments facilitate x-space approximations in the 
style of, e.g, Ref. [35] which are sufficient for most phenomenological purposes, including the 
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determination of the third-order QCD corrections to the Paschos-Wolfenstein relation [36] used 
for the extraction of sin 2 Q\y from neutrino DIS. 

The outline of this article is as follows. In Section[2]we briefly specify our notations and write 

(2) 

down the complete second-order results dca (x) for the above coefficient-function differences. 
We discuss their behaviour at the end points x = and x = 1, and provide compact but accurate 
parametrizations for use in numerical applications. We then proceed, in Section [3] to our new 

(3) (3) 

results for the five lowest odd moments of dc 2 L and even moments of 5c 3 , as a byproduct deriving 
the third-order coefficient-function correction to the Gottfried sum rule. These three-loop moments 
are presented in a numerical form and employed to construct x-space approximations valid at 
x > 10~ 2 . In Section |H we address the numerical implications of our results. In particular we 
discuss the higher-order QCD corrections to the Paschos-Wolfenstein relation. Our findings are 
finally summarized in Section [51 The lengthy full expressions of the new third-order moments in 
terms of fractions and the Riemann ^-function can be found in the Appendix. 



2 The complete second-order results 

We define the even-odd differences of the CC coefficient functions C a for a = 2, 3, L as 

6C 2 , L = CZ£* P -CZ£* P , 5C 3 = C 3 v ^-C 3 Vp+Vp . (2.1) 

The signs are chosen such that the differences are always 'even - odd' in the moments N accessible 
by the OPE [26], and it is understood that the d abc d a b c part of C^ p+Vp [27, 31] is removed before 
the difference is formed. The non-singlet quantities (|2.1I) have an expansion in powers of oc s , 

8C fl = J> s z 5c? (2.2) 

1=2 

where, as throughout this and the next section, we are have normalized the expansion parameter as 
a s = oc s / (4tt). There are no first-order contributions to these differences, hence the sums start at 
I = 2 in Eq. (TO) . 

All known DIS coefficient functions in massless perturbative QCD can be expressed in terms 
of the harmonic polylogarithms H m ,...,m K {x) with mj = 0, ±1. Our notation for these functions 
follows Ref. [37] to which the reader is referred for a detailed discussion. For w < 3 the harmonic 
polylogarithms can be expressed in terms of standard polylogarithms; a complete list can be found 
in Appendix A of Ref. [23]. A Fortran programs for these functions up to weight w = 4 has 
been provided in Ref. [38], with an unpublished extension also covering w = 5. In the remainder 
of this section we employ the short-hand notation 

#0,...,0,±1,0,...,0,±1,...(X) = H ±(m+l),±(n+l),...( x ) C 2 - 3 ) 
m n 

and additionally suppress the arguments of the harmonic polylogarithms for brevity. 
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Exact expressions for (moments of) the coefficient functions will be given in terms of the 
S\J(N C ) colour factors C A = N c and C F = (N 2 - \)/{2N c ), while we use the QCD values C A = 3 
and Cp = 4/3 in numerical results. All our results are presented in the MS scheme for the standard 
choice fi r = = Q of the renormalization and factorization scales. 

(2) (2) 

The second-order coefficient functions 8c 2 and dc L for the even-odd differences of F 2; l read 
5c® (x) = C F [C F -C A /2](-^ + 112(l+x)- 1 ^3 + ^- 1 + ^+^x 2 

n 144 , 

- 4(K 3 + 136^3* + 8^ 2 + 56£ 2 * + 96^ 2 x 2 - —^ 2 x 3 - 32H_ 2;0 

+ 96H_ 20 (1 +X)- 1 + 128H_ 20 x- 128H_ 1 (1 +x)~ 1 ^ 2 + 48H_ 1 ^ 2 

-144H_ 1 ^ + 32H_ li _ li0 -128H_ li _ li0 (l+^)- 1 -224H_ 1; _ 1 ^ 

16 , , 144 . 

+ 64H 1 >0 + — H j >Q x 2 + 64H l fi x + 96H x $x 2 - — H x ^x 3 

- 64H_ li0>0 + 160H_ 1>0|0 (1 + 160H_ 1A0 x+64H_ 1)2 (l 

28 16 292 

- 32H_ ! 2 + 32H_ ! 2* + y H - 32H ( 1 + - 1 - yH^" 1 - _H ;t 

144 

+ 32H (1 +x)-% + — U Q x 2 - 16H ^ 2 + 16H £ 2 * - 16H 00 - 64U x 

9 144 o , 

- 96H 00 * 2 + — Hq o* 3 + 24H 0>0>0 - 48H 0)0)0 (1 + x) 1 - 24H 000 * 

-32H 1 +32H 1 ^-16H 2 -16H 2 ^+16H 3 -32H 3 (1+jc)" 1 -16H 3 ^ , (2.4) 

5c[ 2) W = C F [C F -C A /2](y^ 1 -^ + ^+y^ 2 + 64^ + 32^ + 64^ 2 

96 64 
-y^ 3 +64^^-64^^^- 128H_ l5 _ lj0 ^ + 64H_ l!0 + yH_ l!0 x- 2 

96 32 
-32H_ l!0 x- 1 +64H_ 1;0 ^ + 64H_ l!0 x 2 -— H_ 1;0 ^ 3 +64H_ 1A0 ^+— H 



64 , 448 96 , , 96 ,\ 

• - 1 __H ^+ T H ^ 2 -32H ^-64H ^ 2 + T H ^ 3 J . (2.5) 

.(2) 



The corresponding quantity 8c 3 for the charged-current structure functions F 3 is given by 

&f(x) = 8c< 2) (*) -C F [C F -C A /2](-^ + ^- 1 + ^+^ 2 + 32^ 3 

144 

+ 96^ 3 x - 16^ 2 + 48£ 2 * + 80^ 2 x 2 - — ^ 2 x 3 + 32H_ 2j0 + 96H_ 2j0 x 
-32H_ 1 C 2 -96H_ 1 C 2 x-64H_ 1) _ 1)0 -192H_ lj _ lj0 x + 64H_ li0 

+ yH_ lj0 ^ 2 - leH^!^- 1 +64H_ 1;0 ^ + 80H_ l!0 ^ 2 - ^H_ 1;0 * 3 

+ 32H_ 1)0)0 + 96H_ l Q Q x- yHo^ 1 - ^H - ^pH *+ ^H * 2 

+ 16H 0j0 - 48H 0;0 * - S0H 0fi x 2 + ^H 0;0 x 3 ) . (2.6) 



Expressions equivalent to Eqs. (12.41) and (12.61) have first been published in Refs. [20] and [22], 
respectively, and were later confirmed in Ref. [23]. To the best of our knowledge, on the other 

(2) 

hand, the function 5c L has not been documented in the literature before, see, e.g., Ref. [39] and 
references therein. It was however calculated by the authors of Refs. [20-22], distributed in a 
Fortran package of the two-loop coefficient functions, and employed for the parametrizations 
of Ref. [40]. Our expression (|2.5I) agrees with this unpublished result. 

It is instructive to briefly consider the end-point limits of the above results. Suppressing the 
ubiquitous factor CpCpA = Cp [Cp — Ca/2], the small-* behaviour of Eqs. (12.41) - (12.61) is 

4 In 3 * - 81n 2 Jc-(28-16^ 2 )lnjc - 64+ 8^ 2 +72^ 3 + . . . 
41n 3 x-161n 2 Jc+(12 + 16^ 2 )lnjc + 44 + 24^ 2 + 40^ 3 + ... 

321nx-48 + ... . (2.7) 

Thus the even-odd differences are not suppressed with respect to the vp +vp two-loop non-singlet 
coefficient functions for x — > 0: the same powers of \nx enter Eqs. (12.71) and those quantities. 
At large x, on the other hand, all three functions hc^ 1 are suppressed by factors (1 — x) 2 times 
logarithms, reading 

-(I2-8C2) [1 -x]C F C FA + 0([l-x] 2 ) 
(20- 8 £ 2 ) [1 -x] CfCfa + O ([1 -x] 2 ) 
(32-16^ 2 )[l-*] 2 C F C FA + 0([l-*] 3 ) • (2.8) 

The differences Sc^ (x) and 8c^ (x) (both multiplied by -1 for display purposes) are compared 
to the corresponding even-/V vp+Vp coefficient functions in Fig.CQ The quantities (12.41) and (12.51) 
are negligible at x > 0. 1 and at x > 0.3 , respectively, but indeed comparable to the even-moment 

coefficient functions at small x. The corresponding results for F3 are qualitative similar to those 

(2) 

for F%, but with 8c 3 (x) small down to x ~ 0.01 . 

For certain numerical applications, for instance for use with complex-/V packages like Ref. [41], 
it is convenient to have parametrizations of Eqs. (|2.4I) - (|2.6I) in terms of elementary functions. 
With an error of less than 0.1% these functions can be approximated by 

~ {-9.1587- 57.70^ + 72.29 x 2 - 5. 689x 3 -xL (68. 804 + 24.40Lo 

+ 2.958L 2 , ) + 0.249 L + 8/9 Lq (2 + Lo) }(1—x) , 
~ {-29.65 + 116.05x-71.74x 2 -16.18jc 3 +xL (14.60 + 69.90x 

- 0.378 L 2 Q ) - 8.560L + 8/9 Ll (4 + Lo) } (1 -x) , 
~ { 10.663 - 5. 248 x - 7. 500 x 2 + 0.823 x 3 +xL ( 11. 10 + 2.225 L 

-0.128L 2 ) )+64/9L }(l-Jt) 2 . (2.9) 

Here we have employed the short-hand Lq = \nx and inserted the QCD values of Cp and Ca- 



8c 2 (x) ~ — 
8c^ (x) ~ — 

Sc^ (x) ~ — 



dc [ 2 2) (x) 
bcf\x) 
5c[ 2) W 



dc^ix) 
bcf\x) 
8cf(x) 
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-4 -3 -2-1 -4 -3 -2 

10 10 10 10 „ 10 10 10 10 



Figure 1: The odd — even non-singlet differences —hc 2 L (x) of Eqs. (12.41) and (12.51) . compared at 
x < 0.8 to the corresponding even-Af coefficient functions calculated in Refs. [19, 20, 23]. 



3 Third-order moments and approximations 

Recently the first five odd-integer moments have been computed of the third-order coefficient func- 
tions for F^l Vp in charged-current DIS, together with the corresponding moments N = 2, . . . , 10 

for F^ p ~ yp [28]. Unlike previous fixed-Af calculations, the complete three-loop results for F 2 V £ vp 
[24,250 and F^ p+vP [27] facilitate analytic continuations to these values of N. We have performed 
this continuation using the x-space expressions in terms of harmonic polylogarifhms [37] and the 
Mellin transformation package provided with version 3 of Form [42]. Thus we are in a position 
to derive the respective lowest five moments of the hitherto unknown third-order contributions to 
the even-odd differences (12 .11 ). These moments represent the main new results of this article. With 
one exception (see below) the exact SU(N C ) expressions are however deferred to the Appendix. 

Here we present numerical results for QCD, using the conventions introduced at the beginning 
of Section (2[ recall especially a s = a s /(4n) and the scale choice [i r = = Q- I n addition rij- 
denotes the number of effectively massless quark flavours, and we use the notation dC a ,N for the 
A^-th moment of dC a (x). The results for F2 and Fi read 

8C 2 ,i = -4.378539253^ + ^ (-125.2948456- 0.6502282123 /y) 
5C 2i3 = -0. 138066958 a 2 , + c£ (-5.554493975 + 0.1939792023 

1 The 0Cs coefficient functions for this process are those of photon-exchange DIS, but without the contributions of 
the flu flavour classes, see Fig. 1 of Ref. [25], where the two photons couple to different quark loops. 
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and 



5C 2 , 5 = 


-0.032987989^ 


+ «s(" 


-0.707322026 + 0.0004910378^; 


6C 2 ,7 = 


— 0.013235254a s 


+ «s (" 


A AAO 01£0^ A A1A1 A£ AiC/^A ~ 1 

-0.008816536 — 0.0201069660/^ 


5C2,9 = 


-0.006828983a* 


+ «s( 


0. 133159220 - 0.0200289710/^; 


5C L i = 


-2. 138954096 a 2 


+ a 3 f- 


- 106.6667685 + 3.294301343 n f ) 


5C L , 3 = 


-0.078259985^ 


+ a, (- 


-9.239637919 + 0.2718024935^) 


5C L , 5 = 


-0.016892540a* 


+ (- 

1 S V 


-2.548566852 + 0.0650677 125 n f ) 


5C L , 7 = 


-0.0062631 13 a 2 




- 1 .075400460 + 0.025 1053847 n f ) 


8Q,,9 = 


-0.00300123 la 2 


+ «?(- 


-0.560603262 + 0.0122952192/^) 



(3.1) 



(3.2) 



The lowest even moments for the structure function are given by 

8C 3 , 2 = -0.1135841071a 2 + a 3 s ( 8.386266870 + 0.0605431788/^) 



8C 3 ,4 = -0.0683669250< + < (-1.237248886 + 0.0971522112^) 

8C 3i6 = -0.0350849853 a 2 + a] (-1.370404531+0.0496762716/^) 

8C 3 ,8 = -0.0208455457 a 2 + a 3 (-l .052847874 + 0.0282541 123 n f ) 

8C 3j i = -0.0137316528a 2 + a 3 (-0.798850682 + 0.0177100327/^) . (3.3) 

The new a 3 contributions are rather large if compared to the leading second-order results also 
included in Eqs. (|3.1I) - (|3.3I) with, e.g., a s = 1/50 corresponding to a s ~ 0.25. Except for the 
lowest moment for a = 2,L, on the other hand, the integer-TV differences dC a ^ are entirely negli- 
gible compared to the Vp±Vp moments of Refs. [28,31]. 

Before we turn to the jc-space implications of Eqs. (13.11 ) - (13.31 ), let us briefly discuss some 
interesting structural features of our third-order results. For this purpose we consider the exact 

-2 



(3) 

SU(/V C ) expression for the lowest moment of 8c 2 given by 



(3) _ _ 2 ( 175030 49216 404720 562784 2 33200 
8c 21 - C F C FA ^— _ fe + __Q3___c a + _£ 2fe 

4160 8992 3 1472 2 

2 / 303377 41350. 363896. 396824 2 26000. 

+ Cf CfA \—m~ + ~~rf~ 8l~^ 3 + IsT^ 2 9~^ 3 

+ 25616 1^6 56432 A ^ 



9 " 3 3 ^ 315 

, „.„6 3056. 39592. r r _ 

+ C F C^ ^_ - — C2 + ^t-C 3 + — ^3 " "135-C 



^8786 3056. 39592 „ 1408. . 30424 „ 2 1792, 



(2) 

As all other calculated moments of the functions Sc^ (x), this result contains an overall factor 
CfA = Cp — Ca/2 = —l/(2N c ). Hence the third-order even-odd differences are suppressed in the 
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large-/V c limit as conjectured, to all orders, in Ref. [34] on the basis of two-loop results in particular 
for N = 1 Adler and Gottfried sum rules, for a recent discussion see also Ref. [43]. In fact, up to 
the additional fin contribution absent in charged-current DIS (recall Footnote 1), 

A (3) d abc d ahc ( oog 1472, 256 2 1280 

A e.m.C 2 j = (^-288+96^2 + ^-^-— ^ 2 -^-^ 5 

= -33.67693293^ in QCD , (3.5) 

Eq. (13.41) represents the third-order coefficient-function correction to the Gottfried sum rule (gsr3 
since the Adler sum rule involving the non- singlet coefficient function C2,i of the vp — vp combi- 
nation does not receive any perturbative or non-perturbative corrections, see, e.g., Ref. [44]. 

Another interesting feature of the functions Scj^ 3 in Eq. (|2.2I) is the presence of ^-functions 
up to weight 21 in the integer moments, e.g., terms up to ^| an d occur in the third-order result 
(13.41) . This is in contrast to the 'natural' (OPE-based) moments of C v a p± ^ p which only include 
contributions up to weight 21 — 1, see Refs. [28-31]. Yet the x-space expressions of all these 
quantities consist of harmonic poly logarithms up to weight 2/ — 1 corresponding to harmonic sums 
up to weight 21. Note also that, in the approach of Refs. [20-22], the absence of weight-2/ terms 
in the natural moments appears to require a cancellation between different diagram classes. 

We now return to the numerical moments (|3.1I) - (13.31) and investigate their consequences for 

(3) 

the x-space functions dc a (x). We follow an approach successfully used, for instance, in Ref. [35] 
when only the coefficient-functions moments of Refs . [29-3 1 ] were known. Based on the two-loop 
end-point behaviour in Eqs. (|2.7I) and (12.81) we expect small-.* terms up to In 5 x and In 3 x in Sc^l (x) 

(3) ? ' 

and 5c L (x), respectively, and large-* limits including contributions up to (1 — x)^ a In ( 1 — x) with 

(3) 

r| 2 3 = 1 and r\ L = 2. Thus the x-space expressions of 5c a will be of the form 

fi\ / 2 Inr \ 7 ~ 2r la 

5cl 3) (x) = (1 -jc)* ( £ A m ln m (l -x) + 5c a smooth W + *i ■!=-) + £ B n \n n x (3.6) 

m=l n=l 

where the functions 8cJ mooth (A:) are finite for < x < 1. For moment-based approximations a sim- 
ple ansatz is chosen for these functions, and its free parameters are determined from the available 
moments together with a reasonably balanced subset of the coefficients A m and B n . This ansatz and 
the choice of the non-vanishing end-point parameters are then varied in order to estimate the re- 

(3) 

maining uncertainties of 8c« (x) . Finally for each value of a two (out of about 50) approximations, 
denoted below by A and B, are selected which indicate the widths of the uncertainty bands. 

For F2 and Ft these functions are, with Lq = \nx, x\ = 1 — x and L\ = \i\x\, 

5c^iW = (54.478L 2 + 304.6Li + 691.68x)jti + 179.14L -0.1826L 3 J 

+n / {(20.822x 2 -282.1(l + f))xi -(285.58*+ 112.3- 3.587 L§)L } , 

dcflix) = -(13.378L 2 + 97.60Li + 118.12x)jci-91.196L^-0.46444 (3.7) 
+ fly{(4.522Li +447.88(1 + f)) x\ + (514.02*+ 147.05 + 7. 386Lo)L } 



2 Note that our overall normalization and expansion parameter differ from those of Ref. [34]. Consequently the 
corresponding GSR coefficients (13. U . (13.41 ) and ( 13.5b are larger by a factor 4 z /3 at order oc s ' than in their notation. 
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and 



(3) 



= _ (495.49 x 2 + 906.86)* 2 - 983. 23xx l L + 53.706 L 2 Q + 5. 3059 4 
+ n f {(29.95jc 3 - 59.087* 2 + 379.91) x\ - 273.042xL§ + 71.482jciLo} 



8c[^(x) = (78.306Li+6.3838x)x 1 z + 20.809xciL -l 14.47 Lg-22.222Lg 

+n / {(12.532Li + 141.99* 2 - 250.62*) x\ - (153.586* - 0.6569) jqL } 

The corresponding results for F3 read 



(3.8) 



S4 3 i(x) 



= (3.216L 2 + 44.50L! - 34.588) xi + 98.719L 2 , + 2.6208 4 
-n / {(0.186Li + 61.102(l+x))xi + 122.51xL -10.914L^-2.7484} , 

= -(46.72L 2 + 267.26Li +719.49 x)x\ - 171.98L + 9.470L^ 
+ n f {(0.8489Li + 67.928 (1 + f )) jq +97.922xL - 17.070L 2 , - 3.132L^} 



(3.9) 



(3) 

The resulting approximations for the vp — vp odd-TV coefficient functions c 2 L {x) are compared 
in Fig. [2] to their exact counterparts [24, 25] for the even-TV non-singlet structure functions. The 
third-order even-odd differences remain noticeable to larger values of x than at two loops, e.g., up 
to x ~ 0.3 for F2 and x ~ 0.6 for Fi for the four-flavour case shown in the figure. The moments 
TV = 1, 3, . . . , 9 constrain dc^ L (x) very well at x > 0.1, and approximately down to x ~ 10~ 2 . 

4000 



2000 - 



-8000 



H — I I I I I I I 



1 1 — I I I I I I I 



H 1 — I I II I 



cfi(x,n f =4) 




even N 

odd N (A, B) 



J I I I I I I 1 1 1 1 



500 - 



-500 



-1000 



"T 1 — I I I I I I I 



"T 1 — I I I I I I I 



I 1 — I I I I 



Ox,n, = 4) 




even N 
odd N (A, B) 



-1500 h / 

_£J I I I I I I I I I I I 



10 



10 



10 



10 



10 



10 



X 



Figure 2: The exact third-order coefficient functions of the even-TV structure functions F2][ +Vp 
[24, 25] for four massless flavours, and the corresponding odd-moment quantities obtained from 
these results and the approximations (13.71) and (13.81) for the even - odd differences. 
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For some applications, such as the Paschos-Wolfenstein relation addressed in the next section, 

(3) 

one needs the second moments of the functions dc 2 [(x). These quantities can now be determined 
approximately from the above x- space results, yielding 

8c§2 = -20.19±0.39 + (0.69^0.040)/^ 

5c[ 3 ^ = -24.75 ±0.15 - (0.792 ±0.014)^. (3.10) 

Here the central values are given by the respective averages of the approximations A and B in 
Eqs. (13.71) and (13.81) which directly provide the upper and lower limits. 

Returning to x- space we recall that uncertainty bands as in Fig. [2] do not directly indicate the 
range of applicability of these approximations, since the coefficient functions enter observables 
only via smoothening Mellin convolutions with non-perturbative initial distributions. In Fig. [3] we 
therefore present the convolutions of all six third-order CC coefficient functions with a characteris- 
tic reference distribution. It turns out that the approximations (13.71) and (13.81) of the previous figure 
can be sufficient down to values even below x = 10 . The uncertainty of 5c 3 (x), on the other 
hand, becomes relevant already at larger values, x < 10~ 2 , as the lowest calculated moment of this 
quantity, N = 2, has far less sensitivity to the behaviour at low x. 




j i i i i i i i I i i i i i i i i I i i i i i i i i 



-3 -2 -1 

10 10 10 1 

X 

Figure 3: Convolution of the six third-order CC coefficient functions for ^2,3,/, m Yp ±Vp [24,25, 
27] and \p — vp [Eqs. (|3.7I) - (13.91 )1 DIS with a schematic but typical non-singlet distribution/. 
All results have been normalized to f(x), suppressing a large but trivial variation of the absolute 
convolutions for small and large values of x. 
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The three-loop corrections to the non- singlet structure functions are rather small even well 
below the x- values shown in the figure - recall our small expansion parameter a s : the third-order 
coefficient are smaller by a factor 2.0 • 10~ 3 if the expansion is written in powers of a s . Their sharp 
rise for x — > 1 is understood in terms of soft-gluon effects which can be effectively resummed, if 
required, to next-to-next-to-next-to-leading logarithmic accuracy [45]. Our even-odd differences 

(3) 

8ca (x), on the other hand, are irrelevant at x > 0.1 but have a sizeable impact at smaller x in 
particular on the corrections for Fi and Ft- 



4 Applications 

The approximate results for dc a (x) facilitate a first assessment of the perturbative stability of the 
even-odd differences (|2.1I) . In Fig.|4]we illustrate the known two orders for F2 and Ft for oc s = 0.25 
and «y = 4 massless quark flavours, employing the same reference quark distribution as in Fig. [3] 
Obviously our new a s 3 corrections are important wherever these coefficient-function differences 
are non-negligible. On the other hand, our results confirm that these quantities are very small, and 
thus relevant only when a high accuracy is required. Presently this condition is fulfilled only for 
the determination of the weak mixing angle Qw from neutrino DIS to which we therefore turn now. 




-4 -3 -2 -1 -4 -3 -2 -1 

10 10 10 10 10 10 10 10 

X X 



Figure 4: The first two approximations, denoted by LO and NLO, of the differences (|2.2I) for 
F2 and Ft in charged-current DIS. The results are shown for representative values of a s and ny 
after convolution with the reference distribution /(jc) also employed in Fig. 3. The dashed curves 
correspond to the two approximations in Eqs. (13.71) and (13.81) for the new oc s 3 contributions. 
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For this purpose one considers the so-called Paschos-Wolfenstein relation defined in terms of 
a ratio of neutral-current and charged-current cross sections for neutrino-nucleon DIS [36], 



R 



o(v^ -> fi-X) - o(v^N -> /u+X) 



(4.1) 



R~ directly measures sin 2 Bw if the up and down valence quarks in the target carry equal momenta, 
and if the strange and heavy-quark sea distributions are charge symmetric. At the lowest order of 
perturbative QCD one generally finds 



R 



LO 



1 

2 



sin 2 



w ■ 



(4.2) 



The quantity (14.11) has attracted considerable attention in recent years due to a determination of 
sin Qw by the NuTeV collaboration [11]: within the Standard Model their result is at variance 
with other measurements of this quantity [1], see also Refs. [13-15] for detailed discussions. 

Beyond the leading order Eq. (14.21 ) receives perturbative QCD corrections which involve the 
second moments of coefficient functions for the vN — \N neutral- and charged-current structure 
functions^] Armed with the results of Sections [2] and [3] we are now able to finalize the corre- 
sponding a 2 contribution for massless quarks [14] and to present an accurate numerical result at 
order oc s 3 . We denote by q~ = q — q the second Mellin moments of the valence distributions of the 
flavours q 



q = I dx x{q{x) — q{x)) . 
JO 



(4.3) 



The QCD corrections to R can be expanded in inverse powers of the dominant isoscalar combina- 
tion u + d~ of the parton distributions - recall that the measurements of this ratio are performed 
for (almost) isoscalar targets. After inserting the expansion of the MS coefficient functions in 
powers of oc s , the Paschos-Wolfenstein ratio Eq. (14.11) can be written as 



R- 



7 U 

Sr + ~ 



d +c 



+ (g 2 L-gl) 



u +d 

8 a s a 
— - H '■ 

9 n n- 



3{gLu-gRu) + {gLd-8Rd) 



+ 



5175965 356 



52488 729 



15127 
1944 " 
586 

S3 



27 



89 61 
8P 2 + 27^ 3 
128. 2 190. 

<2+ — £ 5 



32 , 2 
45 ^ 



83 

T62 n / 



405 



81 



9062 2 „ 
729 / 3 /S3 



226 , 

H ni 

729 f 



-5c (3) + 
32 dC2 ' 2 + 



— 6c (3) ' 
128 5Cl - 2 



}) + {("~+d-)- 2 ) + 0(a 4 s ) . (4.4) 



Here the left- and right-handed weak couplings g^ u , gu, gn u and are related to the weak 
mixing angle sin 2 Q w by 



2 2 2 

8L = 8Lu+8Ld 



- — sin 2 Bw + - sin 4 Bw 



2 2 

8r u + §Rd 



5 • 4 fi 

-sin 0^ 



(4.5) 



Specifically the ratio R includes, besides all v/V — v/V CC coefficient functions, the neutral-current quantity C 3 

3 



NC 



which is equal to its charged-current counterpart C^ N vN at the perturbative orders considered here. 
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Beyond the tree level, of course, these relations receive electroweak radiative corrections, see, e.g., 
Ref. [46]. Eq. (14.41) shows the well-known fact that the relation (|4.2I) receives corrections if the 
parton content of the target includes an isotriplet component, uT 7^ d~, or a quark sea with a C-odd 
component, s~ 7^ or c~ 7^ 0. Notice also that perturbative QCD only affects these corrections. 

The exact second-order contribution in Eq. (14.41) differs from the result in Ref. [14] where the 
function bc L (x) of Eq. (12.51) was not included. The third-order corrections can now be completed 
in a numerical form, using our approximations (13.101) for the second moments of §c 2L {x). For 
rijr = 4 flavours (and disregarding electroweak corrections) we obtain 

R = --sin 2 w + u - + d - U-^sm 2 Q w + f - -sin 2 w J - 

| — [1 + 1.689 Os + (3.661 ±0.002) a 2 ] 1 + o ((u~ + d~)~ 2 ) + o{a 4 s ) . (4.6) 

The perturbation series in the square brackets appears reasonably well convergent for relevant 
values of the strong coupling constant, with the known terms reading, e.g., 1 + 0.42 + 0.23 for 
0C S = 0.25. Thus the a 2 and 0^ contributions correct the NLO estimate by 65% in this case. On 
the other hand, due to the small prefactor of this expansion, the new third-order term increases the 
complete curved bracket in Eq. (14.51) by only about 1%, which can therefore by considered as the 
new uncertainty of this quantity due to the truncation of the perturbative expansion. Consequently 
previous NLO estimates of the effect of, for instance, the (presumably mainly non-perturbative, see 
Refs. [47-49]) charge asymmetry of the strange sea remain practically unaffected by higher-order 
corrections to the coefficient functions. 



5 Summary 

In this article we have presented new results for the coefficient functions of inclusive charged- 

current DIS in the framework of massless perturbative QCD. We have filled a gap in the two-loop 

(2) 

literature by writing down the corresponding difference dc L (x) of the Vp+Vp and vp—Vp 
structure functions Ft . Our main results are the lowest five (even- or odd-integer) Mellin moments 

(3) 

of the third-order corrections §c a (x) for all three structure functions F a= 2 3 l and approximations 
in Bjorken-x space based on these moments which are applicable down to at least x < 10~ 2 . As a 
byproduct we have calculated the related third-order coefficient-function correction to the Gottfried 
sum rule in photon-exchange DIS. 

All our third-order results are proportional to the 'non-planar' colour factor Ca — 2Cf, thus 
confirming a conjecture by Broadhurst, Kataev and Maxwell on the 1/A^ 2 suppression of these 
coefficient-function differences in the limit of a large number of colours A^. Numerically our 
corrections prove relevant in particular for F2 and Ft wherever the differences of the vp + vp and 
vp — vp coefficient functions are not negligible. We have employed the above results to derive the 
second- and third-order QCD corrections to the Paschos-Wolfenstein ratio R used to determine 
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the weak mixing angle from neutrino deep-inelastic scattering. The uncertainty due to uncalculated 
higher-order coefficient functions has been reduced to a level amply sufficient for the foreseeable 
future, i.e., 1% for the coefficient-function factor multiplying the quark-distribution asymmetries. 

Form files and Fortran subroutines with our results can be obtained from the preprint server 
|http: //arXiv . org] by downloading the source of this article. Furthermore they are available from 
the authors upon request. 



Note added 

(3) (3) 

While this article was finalized, the 1 1-th moments of the functions 8c 2 (x) and 5c L (x) have been 
computed [50]. Both results fall into the bands generated by the respective jc-space approximations 
in Section [3l thus confirming the reliability of these uncertainty estimates. 
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Appendix 



(3) 

Here we present the analytic expressions for the Mellin-space coefficient-function differences 5c a N 



which were given numerically in Eqs. (13.11) - (13.31) . We use the notations and conventions as 
specified at the beginning of Section [2] and above Eq. (12.41 ). 

The first moment of Sc 2 (x) has been written down in Eq. (13.41 ) above. The remaining known 
moments of this quantity are given by 

(3) 2 / 1805677051 2648 10093427 1472 2 7787113 

^2,3 - C F C FA ^ 466560 g~<* + 810 ^ 3 3"* 1944"^ 

55336 378838 2 8992 3 

2 / 5165481803 40648 9321697 1456 2 8046059 

+ CfCfA V 1399680 + ^ 5 ~ "SIO - ^ 3 + "l"^ 3 + ~l944~^ 2 

jno ,, , 798328 r2 56432 r3 \ 
- 4984^3 + ^r8-^T$ ) 

/20396669 1792 405586 139573 1408 50392 2 \ 

+ n f C F CpA 1^16640 9 _ ^ 5 + ^05 _ ^ 3 486 _ ^ 2 + "9 _ ^ 3 "^35~^ 2 J ' 
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/ 18473631996593 17584 149815672 1472 2 

F FA V 3827250000 45" 7875 ^ 3 3~~ ^ 3 

291199027 330416 2577928 2 8992 3 

50625 ^ 2 + 45 ^ 3 225 ^ 2 ~ ~&f^ 2 

2 / 16016244428419 47560 1270840912 1456 2 

+ F FA \ 3827250000 + 9 ^ 5 70875 ^ 3 + 3 ^ 3 

1321405949 89128 26658224 2 56432 3 

+ 202500 ^ 2 15 - ^ 2 ^ 3 + 3375 ^ 2 ~ 315 ^ 2 

/ 181199822513 1792 6514448 1652773 
+ n f C F C FA ^ 765450000 9-^5 + ^25~^3 " -3375-fa 

1408 11888 2 
^~ 27" 

2 / 177036089007294328733 27248 65397081433 
F FA V 32934190464000000 63~ ^ 5 + 2646000 ^ 3 
1472 2 340303364748629 2563996 
3~ ^ 3 46675440000 ^ 2 + 315 ^ 3 



4570738447 2 _ 8992 3 
330750 ^ 2 ~ ~63~^ 2 
2 C / 213694072871074531 1821772 438487320707 
+ F FA \ 45177216000000 + 315 ^ 5 18522000 ^ 3 
1456 2 418808510000479 2071492 

+ ~T~ ^ + 46675440000 ^ 2 315 ^ 3 

6241478743 2 56432 3 
+ 661500 ^ 2 ~ 315 ^ 2 

/ 38079608000704561 1792 22115039 
+ n f F FA \ 117622108800000 9~ 13230 ^ 3 

113587875043 1408 2296328 2 

166698000 ^ 2 + ~9~^ 3 4725 ^ 2 

2 / 5676515460744370321603 25664 11165079556403 
F FA \ 1000376035344000000 63~ 375070500 ^ 3 
1472 2 8178803099431493 1648352 

3~ ^ 3 945177660000 ^ 2 + 189 ^ 3 
23488033336 2 8992 3 
1488375 ^ 2 ~ 63 ^ 2 
2 / 32102287673972370020989 1162796 89153747611 
+ F FA \ 6002256212064000000 + 189 ^ 5 3087000 ■ 
1456 2 342078312478997 1332820 
+ ~T~ ^ 3 + 30005640000 ^ 2 189 ^ 3 

3187232017 2 56432 3 
+ 297675 ^ 2 ~~ 315 ^ 2 
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/2 1832 132 134852204299 1792 .... 

+ n f C F C FA ^ 52400649470400000 9~^ 5 + " 1 ^ c ^ 



6271692134 
3274425 ^ 3 
193! 8 24 297943,.., 1408. ... 164116, 



2250423000 ^ 2 + 9 ^ 3 315 ^ z ) ' (A ' 4) 
The corresponding lowest five odd-integer moments for the longitudinal structure function read 

+ Cfe (- ™» - !f C + + ^ + 288^3 + ^f) 

„ „ /1366 496 „ 328 r 224 „ 2 \ 
+ «/C F C FA ^- -^-—^-—gj , (A.5) 

( 3) 2 / 12350749 52516 47 7544 2 

54 j = C F C FA (^ ___ + 352C 5 + 27 C2 + 96^3-^-^2 

2 / 10152961 16412 1168 2 

+ C F C FA I 1296Q 368Q 5 — C, 3 - 242C, 2 + 144^^ + —^^2 

„ „ ( 16757 2936. 16 . 368. 2 \ 

(3) 2 f 735306721 1888 558244 442783 

*x,s " QC^- 17010000 r^ + ^Ts - ^ 3 675"^ 

+ 448^"^?) 

2 ^ 6741265367 736 1285168 51493 

+ CfCfA \ 10206000 ~ " _ 945 _ ^ 3 + ~405~^ 2 

_ . 69608 . 2 \ 
+ 96^3 + ^5-^ 2 J 

^ ^ / 2107157 8816. 11992. 736. 2 \ 
+ "^ C H~ ^5150 + -^{,-—5,-—^ , (A.7) 

,3, = a /354522585410107 _ 47266403 _ 1095179473 

' F FA \ 666792000000 S5 23625 S3 945000 S2 

2 ( 11388456807174161 176925641 4569363329 

+ F FA \ 28005264000000 ^ 5 132300 ^ 3+ 13230000 ^ 2 

y 663878. 2 \ 

+ 72 ^ 3 + ^625-^J 

Z369546282989 124282 220747 184 2 \ 
+ n f C F C FA ^ 50009400000 + - ^250-^2 " ^2 ) , (A-8) 
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§J 3 ) = C C 2 [ — — — " ' ' C I : "" r 

L ' 9 F \ 1323248724000000 5 S5 6125 S3 

325373958827 „ 5184„ „ 296736, 



/ 134645491 1003496947 10528 „ 13247918. 



208372500 ^ 2 + _ 5~~ ^ 3 875 ^ 
2 / 17693872049573089 736 125991917 
+ F FA \ 73513818000000 5~^ 5 99225 ^ 3 
10496201057 288 1688888 2 

+ 23152500 ^ 2 + ~5~^ 3 + 7875 ^ 2 
/ 23852323249607 1 249264 1542176 736 2 \ 
+n f C F C FA \^ 1444021425000 +T7^^ 3 "^3075~ fe "^75"^J ' (A>9) 

Finally the analytic expressions for the moments dc^ N in Eq. (13.31) are 
*(3) ^2/ 840949 9344 650360 1472 2 712328 

47920 30416 2 8992 3 N 
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I 30416 2 8992 3 \ 

-^3 + —^ + -^-^ 



2 / 15979879 38416 580504 1456 2 742390 
+ C f Cfa \—&H + ~^T~ r ^ 243~^ 2 

„ 577264^ 56432 r3 \ 
•4368C 2 C 3 —^-^J (A.10) 

„ „ ( 119522 1792,. 57128,. 46112^ 1408,. Y 40024 2 

+ n f C F C/74 1 L; Co CoCtH C? 

f f pa ^ 729 ^ 9 S5 81 S3 243 S2 9 S2S3 135 S2 

(3) 2 f 21230721185377 23704 292322783 1472 2 

6c 34 - C F C FA (^ 4374000000 + ^5~^ 5 " 20250 ^ + —^3 

5644168873 100792 6477802 2 8992 3 

+ 1215000 ^ ~ ^ 3 + ~&75~ ^ + "6T^ 2 

/ 19991706724601 _ 272933467 _ 1456 

F * A \ 4374000000 S5 20250 S3 3 S3 

6307524619,. 253064,. , 2500616 yl 56432 ,, 3 
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1215000 " z 45 375 ^ 315 

15339664501 1792 755894 21942049 
72900000 + "9" ^ 5 ~ 675 ^ 3 + 60750 
1408,. „ 53128 



^3 + ^^), (A.11) 



9 135 



(3) _ r r 2 ( 172761364527374293 21200 3380925064 1472 ; 
C3 - 6 ~ F FA V 32162295375000 + 63 ^ 5 165375 ^ 3+ 3 ^ 3 
_ 147865501939 2395856 75351016 2 8992 3 

+ 24310125 ^ 2 315 ^ 3 + 6125 ^ 2 + 63 ^ 2 
2 /313157547783370669 1810712 67828543996 1456 2 
+ C f C F a I 64324590750000 315 ^ 5 + 3472875 ^ 3 3~ ^ 3 
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3604225183081 667864 1397140016 2 56432 

486202500 ^ 2 + 105 ^ 3 165375 ^ 2 + 315 ^ 

/ 503591542653161 1792 1 6004944 23420609 
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!„ „ 2135888^2 \ 
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2 / 126830527574348410837327 5795836 4175977929883 
+ F FA V 24009024848256000000 945 ^ 5 + 166698000 ■ 

1456 2 194854342276579 6509876 

3~^ 3 20003760000 ^ 2 + 945 ^ 3 
58805551031 2 56432 3 
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2 ( 2924815993615556996346598663 210944 
F FA V 483334682604559632000000 + 385 ^ 5 

3080312718428437 1472 2 13720175530646448109 



99843767100 ^ 3 ^ 1537594013340000 
8455904 3568998808 2 8992 3 



945 218295 ^ 63 

2 r ^61916581373996975119251441821 66873844 
+ F FA \ 10633363017300311904000000 10395 ^ 5 
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+ 998437671000 ^ 3 3~^ 3 263587545144000 ^ 

75900812 1995571648453 2 56432 3 
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